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Abstract 

Variational Monte Carlo calculations have been made for the ^)(0 hypernu¬ 
cleus using realistic two- and three-baryon interactions. A two pion exchange 
potential with spin- and space-exchange components is used for the AN poten¬ 
tial. Three-body two-pion exchange and strongly repulsive dispersive ANN 
interactions are also included. The trial wave function is constructed from 
pair- and triplet-correlation operators acting on a single particle determinant. 
These operators consist of central, spin, isospin, tensor and three- baryon po¬ 
tential components. A cluster Monte Carlo method is developed for noncentral 
correlations and is used with up to four-baryon clusters in our calculations. 

The three-baryon ANN force is discussed. 
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I. INTRODUCTION 


In this paper we initiate a variational Monte Carlo study of hypernuclei using realistic 
two- and three-baryon interactions involving a A and nucleons. In the past, variational cal¬ 
culations of the s-shell hypernuclei, |]^,0 a few p-shell hypernuclei using appropriate models 
1^ and A binding to nuclear matter were performed using mostly simplihed and cen¬ 
tral nucleon-nucleon (NN) interactions. The aims of these calculations have largely been to 
deduce information about the A-nucleon (AN) and A-nucleon-nucleon (ANN) interactions. 
In addition, as a result of these studies, one could also explore the structure of hypernuclei. 

The reason for using simplihed NN interactions is the hope that the uncertainties in the 
NN interaction will largely cancel out in these calculations. This is because the A-separation 
energies 5 a, which are the main experimental input in these calculations, are the differences 
of the energies of hypernuclei and their cores, i.e., - B\ = ~ where ^E is the total 

energy of the hypernucleus and ^~^E is the ground state energy of the core nucleus. However, 
interactions generate strong NN correlations in the nuclear wave function. A realistic NN 
interaction will generate central, spin, spin-isospin, tensor and other two-nucleon correlations 
ii- In addition, there are signihcant three-nucleon correlations. In a hypernucleus, because 
of the operator dependence, these correlations may interact in a complicated manner with 
the AN and ANN correlations. This whole group of correlations then interact with the 
two- and three-body operators of the two- and three-baryon interactions. This can result 
in important contributions to the A-binding energy compared to the use of only central NN 
correlations generated by purely central NN interactions. 

In this study, we use realistic two- and three-nucleon interactions and wave functions 
to see their effects on hypernuclei. There is another important aspect of the present study, 
the development of a methodology and variational program for hypernuclei. For the p- 
shell hypernuclei, we adopt the cluster Monte Carlo(CMC) technique developed in Ref. 0 
which we refer to as PWP. As a hrst step, we study ^((O. Generalization to ^aO 'wiH be 
straightforward. Further development of the one-body part of the nuclear wave functions 
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will be needed for other p-shell hypernuclei. We intend to cover a wide range of hypernuclei 
in order to have reliable information on three-baryon forces, because in this study we hnd 
that the role of the three-body ANN interaction is greatly altered from that found in some 
previous studies. 

There are a few calculations of hypernuclei in which realistic NN forces have been used. 
One such calculation is by J. Carlson in which he explicitly considers the AN —EN 
channel in ^He and ^He using the Nijmegen interaction. This study shows that the Nijmegen 
interaction underbinds the four-body hypernuclei and that the five-body hypernucleus is 
unbound relative to a separated a and A particle. Also, it does not reproduce the spin¬ 
splitting in the four-body hypernuclei. To resolve the classical over-binding |]I|,|9|,pT| problem 
of ^He, Bando and Shimodaya 0 , and Shinmura et al. m have also performed calculations 
with Reid soft core and Hamada-Johnston NN potentials by calculating effective interactions 
using a G-matrix approach. All these calculations of the A-seperation energies B\ are 
conhned to s-shell hypernuclei. 

Since no experimental data for exists, we generate ’’pseudo-experimental” or semi- 
empirical data for this system. This is done in Sec. II where we also briefly discuss the 
experimental status of A-seperation energy data. In Sec. Ill we review the A and nucleon 
two- and three-body potentials. Section IV deals with the variational wave function. In Sec. 
V we describe the techniques of our calculations for In Sec. VI we present our results, 
and Sec. VII contains our conclusions. 


II. A-SEPARATION ENERGIES 


If we combine the results of previous experiments summarized in Ref. [^, the in-flight 
reaction (K“, tt”), 00 and the associated production (7r+,K+), we have now almost 
30 well established hypernuclei with a wide range of baryon numbers A < 81 and orbital an¬ 
gular momentum ijy < 3. The hypernuclei that are relevant for the empirical determination 
of the Ra value of are i|C, ^^Si, ^oca, 5|V, and 
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We use here three approaches for the empirical B\ value of In the first approach 
microscopic calculations of i?A of the above hypernuclei were carried out with phenomeno¬ 
logical two- and three-body AN and ANN interactions that were previously obtained from 
studies of Ap scattering, the s-shell hypernuclei, ^Be as a representative of p-shell hyper¬ 
nuclei in the 2a + K model, and the A-binding to nuclear matter |^|-Q . The A-seperation 
energies B\ are obtained from a Schroedinger equation with a A-nucleus potential and 
an effective mass m\ which are obtained in the local-density approximation using the Fermi 
hypernetted chain technique for the A binding to nuclear matter. Such a procedure gives 
a good account of the i?A data of the above hypernuclei for £a < 3. Using this procedure, 
we calculate the difference, ASa, of the A-separation energies of and This gives 
A5 a = 0.30 MeV. 

In the second approach, we use the purely phenomenological technique adopted by Mil- 


lener et al. [jT^. Here we use a Woods-Saxon A-nucleus potential whose parameters were 
fitted to the B\ data of the above mentioned hypernuclei, 

UoA 


UaW = 


l + exp{^) ’ 


( 2 . 1 ) 


where, Vqa = —28.0 MeV, c=(l.128-1-0.439A a=0.6 fm give a very good account 

of the Ha data. This procedure gives AHa=0.40 MeV. 

In the third approach, we consider a density-dependent A-nucleus potential [pT| -pI| of 
the form 


UA{r) = Ap{r) + Bp‘^^^{r), 


( 2 . 2 ) 


where p(r) represents the nucleon density. These are taken from Ref. [^. We choose s- 
and p-orbits in and the p- and f-orbits in to fix the parameters A and B. The 
experimental binding energies are 12.5±0.35 and 2.5±0.5 MeV in aO and 16.0T1.0 and 
2.5±1.0 MeV in Fitting these energies results in a reasonable fit to all binding energies 
from to ®aY for < 3. The resulting AHa is 0.76 MeV. 
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If we combine the results of the above three approaches and at the same time bear in 
mind that the experimental uncertainty in the value of is ±0.35 MeV, we may 
reliably £x the empirical B\ value of as 13.0 ± 0.4 MeV. We shall make use of this value 
in our calculations. Other approaches, such as relativistic mean-held theories or the local- 
density approximation using a Skyrme interaction can also be employed in the empirical 
determination of the B\ value of But, since these approaches are consistent with the 
approaches that we have adopted above, we do not feel that their inclusion will affect our 
results. 


III. HAMILTONIAN 


In an A-baryon hypernucleus, we will consider the hrst A-1 baryons to be nucleons. We 
will use 4/ to refer to the full wavefunction of the hypernucleus and 4 'at to refer to the 
ground-state wave function of the A-1 nucleons. The full Hamiltonian H can be written as 


H = H^ + Ha, 


(3.1) 


where Hm is the nucleon Hamiltonian: 


A-1 ^2 A-1 

fcl 2m 


A-1 

E Vm., 

i<j<k 


and 


(3.2) 


^2 A-1 A-1 

= “oEW^A + X! '^*A +Y1 ^*iA , 

ZrriA 

is the part of the full Hamiltonian due to the A-particle. 

The A-separation energy, B^, of a hypernucleus is then given by 


(3.3) 


< 4/Ar TTat Tat > 

< > 

< 4/ AfI ^/ai > 

< T|T > 


(3.4) 


Our goal is to calculate Ba using a variational principle for the two components of Eq.(3.4). 
In this section we briehy describe the two- and three-body baryon interactions that we have 
used in this study. 
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A. AN potential 


Two-pion-exchange(TPE) is a dominant part of the AN potential, that in turn is mainly 
determined by the strong tensor one-pion-exchange(OPE) component acting twice. More¬ 
over, there is the K-exchange interaction that primarily contributes to the AN exchange 
potential. The tensor part of the AN interaction is very weak because the shorter range K- 
and iP*-exchanges that are responsible for this are of opposite sign and nearly cancel each 
other. (In the case of the NN interaction the vr-exchange and p-exchange tensor components 
do not cancel so completely, because their masses are quite different.) 

We use an Urbana-type [^] potential with spin- and space-exchange components and a 
TPE tail which is consistent with Ap scattering below the S threshold, 

1 


VAN{r) = Vo(^^)(l - e + (^Px) + ^v^T^{r)fTA ■ (Tn, 


(3.5) 


= Vc{r) - vTl{r) , 


(3.6) 


Wr 


Vrir = 


1 + 


where T 7 r(r) is the OPE tensor potential 


j’.=(i+-+4)-(i 


2\2 


(3.7) 


(3.8) 


x=/ir, /i=0.7 fm“^ is the pion mass, and the cut-off parameter c=2.0 fm“^. is the space 
exchange operator and e is the corresponding exchange parameter. The v = (us-|-3ut)/4 and 
= Vg — Vt are respectively the spin-average and spin-dependent strengths, where Vg and Vt 
denote the singlet and triplet state depths, respectively. (Note that following the convention 
of Ref. P|, the Hamiltonian effectively contains +Vc, -v, -t-Uo-, -Vg, and -u*.) Finally, Uc(r) is 
a short range Saxon-Wood repulsive potential. The various parameters are 


Us = 6.33MeV, Ut = 6.1MeV, e = 0.3, lTc = 2137MeV, R = 0.5fm, a = 0.2fm. (3.9) 
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These parameters are consistent with the low energy Ap scattering data that essentially 
determine the spin-average potential v. The parameter e for the space-exchange strength is 
fairly well determined from the A single-particle scattering data For a detailed acconnt 
of the determination of the other parameters see Ref. |^. We point ont that becanse of 
the non-central AN and ANN correlations introdnced in the next section, the AN spin-spin 
potential will have a non-zero contribntion even in a closed-shell system snch as 


B. ANN potential 


When a AN potential that hts the Ap scattering is nsed, the Ra for hypernnclei with A>5 
are almost a factor of two too large. This is an old resnlt that has been confirmed by varions 
analyses. In the present work we also hnd that the nse of a realistic NN interaction does not 
alleviate this over-binding problem. As in the previons stndies, to resolve the over-binding 
problem, we incorporate a three-body ANN interaction. 

We consider here two types of ANN potentials that arise from projecting ont S,A, etc. 
degrees of freedom from a conpled channel formalism. These are the dispersive and the TPE 
ANN potentials designated as and respectively (see Fig. 1). is expected 

to be repnlsive, and, following Ref. [^, we assnme the phenomenological form 

V^mm = ( 3 . 10 ) 

where Wq is the strength of the potential and TT,{ri\) is given by Eq.(3.8). consists 

of two parts corresponding to p and s wave ttA interactions 

VlJjN = Wp + Ws, (3.11) 


where 


M'p = -(y) (^1.T2){A'ia,X2a}, 


( 3 . 12 ) 


Cs(ti ■ T2){(Ti ■ riA)(cr2 ■ r2A)(/rriA -F l)(/xr2A -F l)y(riA)y(r2A) 

(a^ia?^2a)^ 


(3.13) 
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(3.14) 


Y(r) = — (l - e-"-’) , 

{A,B} = AB + BA, (3.15) 

and 

XiA = {(Ti ■ (Ta)Y (riA) + SiAT^{riA). (3.16) 

Here XiA is the one-pion exchange operator, and SiA is the tensor operator. The component 
Ws is qnite weak and, as in previons studies, we neglect it here; we feel that its effect should 
be studied in future work. 

There are theoretical as well as phenomenological estimates for Cp, but for Wq the 
estimates are purely phenomenological. For example, for Wq ~ 0.02, the reduction in the 
A-binding to nuclear matter(using central correlations) is approximately in accord with the 
suppression obtained in coupled channel (AA^ —> J^N) reaction matrix calculations. For Cp 
theoretical estimates give 1~2 MeV; however, the phenomenological values may not lie in 
this region as the results depend sensitively on the cutoff parameter c that appears in Eq. 
(3.8). In the present study, we have obtained results as a function of the values of these 
parameters. 


C. Two- and three-nucleon potentials 


For the nuclear part of the Hamiltonian, we use NN and NNN potentials that have been 
previously used to study various nuclei, including 0. The NN potential contains the 
Erst six terms of the Argonne ni 4 potential and a Coulomb term: 


+ vcoui{rij)0. 

p=i 


Coul 

ij 


(3.17) 


where the operators are 

Oij j ■) ^ j •! ^ j) i ' ) 7 Sij ^ Sij (^l~i ’ 1 ~ , 

Og“' = l(l + T,,i)(l + r.,,). 


(3.18) 





We shall also refer to these operators by the abbreviations c, r, a, ar, t, and tr. In PWP 
it was fonnd that the 7 < p < 14 terms of Argonne ni 4 and the corresponding p = 7,8 
correlation operators gave a net contribntion of only -0.45 MeV/nncleon. We assnme that 
the presence of a A will not signihcantly modify this resnlt and hence we can safely omit all 
potential and correlation operators for p > 7 when compnting ao). 

The NNN potential is of the Urbana type, which consists of dispersive and two-pion- 
exchange terms: 

V.it = + V,%, (3.19) 




(3.20) 


eye 


eye 


^ijk ~ ^2 ^0 ( {^iji ^jk}{'^i-'^ji . [^iji ^jk] '’’j-'^'k] ) , 


(3.21) 


where the sqnare brackets represent the commntator [A,B]=AB-BA. The constants Aq and 
Uq have the valnes -0.0333 and 0.0038 in Urbana model Vll which we nse here. 


IV. THE VARIATIONAL WAVE FUNCTIONS 


We assnme that a good variational wavefnnetion for a hypernnclens with a closed-shell 
nnclear core and a A-particle can be written as 


IT >= 


n(i T UijA T Uijk') 


IT 


'A-1 


n (1+ u^A) 


i=l 


A-1 

^ n (1+ 

i<j 


0 2 


I^J >, 


(4.1) 


i^j >= n ftAh) n /c(ui)0A(rA)Ai<i)^ ^ > 

2=1 i<j 


(4.2) 


Here Uijk represents a three-baryon correlation operator that has the same structure as Vijk^ 


Uijk ^Uijk') 


(4.3) 


UijA ^A^ijA 


(4.4) 
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The Vijk differs from Vijk through the range c of the cutoff functions of Yt,(y) and T^(r). 
The parameter 6 is referred to as e in PWP; we use 6 here to avoid confusion with the e of 
the space-exchange potential in Eq. (3.5). The parameters 6, 6\, c and ca are determined 
variationally. The label IT stands for independent triplet product of 1+Uijk- Thus, 

+ Uijk + UijA) = 1 + '^UijA+ X! + X! 

IT i<j ^<j i<j<k 

i' <.j' <,k' ^i,j 

+ Y. + ■ ■ ■ ( 4 . 5 ) 

i<j<k 

i <.j <fc ^i,j,k 

The neglected terms are of the type This restriction makes the three-body corre¬ 

lations much simpler to use. As is discussed in PWP, the Uijk and UijA should ideally act last 
as in Eq. (4.1). However, this requires considerably more computer time. The improvement 
in the energy of obtained by this was found to be only -0.19(7) MeV/nucleon. In the 
present work we ignore this correction and compute the energies of both and with 
the three-baryon correlations acting Erst. 

Each operator in the two-baryon interaction can induce the corresponding correlation. 
The and the fc are central correlations that are primarily generated by the repulsive 

cores in the two-baryon interactions. For Uij and Ua, we make the following choice 

n 

Uij = Pp'^piUj)Ofj, (4.6) 

p=2 

and 

m 

f4iA = i:«?(riA)OJ,. (4.7) 

The notation STl in Eq- (4.1) represents a symmetrized-product of the non-commuting 
operators UijUjk ■ ■ ■■ Previous studies [^,0 on few-body nuclei and demonstrate that it 
is probably sufficient to use 2< p < 6 in Eq. (4.6). The pair correlation functions fc and Up 
are generated by minimizing the two-body cluster energy using a quenched potential: 

Vij = (4.8) 

p=2 
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The two-body cluster contribution has been minimized for inhnite nuclear matter at Fermi 
momentum kp, with the boundary conditions fdr > d) = 1, and Up{r- > d) = 0, for p = T,a, 
and ar, and Up{r > dt) = 0, for p = t and tr with their first derivatives zero at r = d or dt. 
For the UiA we consider 

UiA = u^{riA)crA ■ cTi + up^{riA)Px- (4.9) 


In the present study, we have omitted the second, i.e., the exchange correlation term in 
Eq.(4.9). Inclusion of this term increases the computation effort by several fold and prelim¬ 
inary results indicate that it gives a small contribution. This will be the subject of a future 
publication. The spin correlation is 


f A _ f A 

^A J t J s 




(4.10) 


where and //'■ are the solutions of Schroedinger equations with quenched AN potentials 
in singlet and triplet states respectively: 


2fJ,AN 




(4.11) 


The potentials Vs{t) are quenched in the two-pion and spin-exchange parts of the central and 
spin channels: 


Vs{r) 

= Vc{r) - (a27r + I aaVa)Tl{r) , 

(4.12) 

vt{r) 

= Uc(r) - (asvr V a„v„)Tl{r). 

(4.13) 


The spin-averaged correlation function is given by: 


J C A 


(4.14) 


The and have been obtained by minimizing the two-body cluster energy for the 

A-binding to nuclear matter with the asymptotic condition f^{r > dA) = 1. 

The (j)A represents a bound-state wavefunction of a A-particle of mass mA moving in a 
Woods-Saxon potential that is bound to a nucleus of mass (A-l)m: 
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l"A(rA) = 


1"a 


(4.15) 


1 + 

The parameters Va, i?A and a a are determined variationally. The Slater determinant 
A consists of orbitals of nucleons of mass m bound to a hypernucleus of mass (A- 

2)m+mA moving in a Woods-Saxon wine-bottle potential 

1 


V(r) = Vs 


1 4 - Q{r-Rs)/as) 




(4.16) 


with Vs, Rs,(is,C(s, and ps as variational parameters. The coordinates of all the one-body 
orbitals are measured from the center of mass of the whole system, thus making Tj and T 
translationally invariant, i.e.. 


= r • — R 

(4.17) 

+ tuaFa 
m{A — 1) -l- rriA 

(4.18) 


V. THE CLUSTER EXPANSION 


We briefly outline the general framework for the cluster expansion of PWP to calculate 
the expectation values of various operators. These expectation values are needed in the 
evaluation of the energy using the variational wave function (4.1). We demonstrate the 
cluster expansion for the two-body NN and AN potentials: 

(d'lEuj + EuA|d') _ N 

(T|T) D ' ^ ’ 

The N and D can be expanded as a sum of n-body contributions 

A-l A-l A-1 A-1 

N ^ ) Tlij -|- ^ ) ^iA T ^ ) ^ijk T ^ ) ^ijA T ‘ ‘ , (5.2) 

i<j i=l i<j<k i<j 


A-l A-l 

D = 1 + dij + diA 

i<j i=l 


(5.3) 
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The expressions for the purely nuclear riij, dij, riijk, etc., may be found in PWP. The 


contributions of clusters containing a A are similar, e.g. 


TT'iA — ^(1 + 'yjA(l + UiA)'j , 


(5.4) 


di \ — <^(1 + (1 + — 1 , 


(5.5) 


nijA — ( (1 + + Ulj)P (1 + 17/a)(1 + ^ja)] bo + + vja] 

xP [(1 + L(jA)(l + Uja)] (1 + Uij){l + UijA) ) — riij — riiA — tija , (5-6) 

where P and P' are permutation operators. In these expressions 

where Tj/ denotes the Tj of Eq. (4.2) without the anti-symmetrization operator. 

The expansions (5.2) and (5.3) for N and D are divergent |]^. We obtain a convergent 
linked-cluster expansion by expressing 


Cij + CiA - 

^ i<j i=l 

The various c,,- and c^a etc. are obtained from the equation 


N = 


A-l A-1 

pj + X! + 

i<j i=l 




by equating terms with the same zj, zA, ijk, ijA, etc. Thus, 


( 6 . 8 ) 


(5.9) 


Qa — 


‘^iA 


1 + diA 


and 


(5.10) 


^ _ ^ijA ^ijidjA T diA') CjA(dp T djA) CjA(dij T diA) T T C^jA)dijA 

1 + dij + diA + djA + dijA 

In the present work, we have used the CEA expansion of PWP so that all clusters of a 
given spin, isospin, and A content are averaged together. 
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VI. RESULTS 


A. Variational parameters 

We made detailed variational parameter searches for two cases: 1) with no ANN potential 
and hence no correlation, and 2) using a ANN potential with Cp = 0.7 MeV and Wq 
= 0.015 MeV. The rest of the Hamiltonian was as described in Sec. Ill for both cases; in 
particular the NNN potential and Uijk correlation were used in both cases. 

For the case with no ANN potential, we found that the optimal values of all of the 
nucleon correlation parameters are the same as was found in PWP for except that 
the well depth, Vg, of the Woods-Saxon potential changes from -49.1 MeV to -48.9 MeV 
to maintain the same p-wave separation energy, 14.0 MeV, with the 17-body instead of 16- 
body reduced mass. The reader is referred to PWP for these parameter values. The optimal 
parameters for correlation terms involving the A are: for the A Woods-Saxon well, 

Va = —28.3 MeV; Ra = 3.2 fm; oa = 0.5 fm, (6.1) 

which gives an s-wave separation energy of 15.0 MeV; for the Ui\ 

C(2n = «o- = 1-0 ; c^A = 2.8 fm. (6.2) 

In the presence of the ANN potential with Cp = 0.7 MeV, Wq = 0.015 MeV, we found 
that only two of the above optimal values had to be changed. These are the quenching 
parameters a in the NN correlation and a 2 n in the NA correlation. The variational energy 
is sensitive to a and we made several searches at other values of Cp to determine 

a = OM - 0.1 Cp, 0 < (Pp < 1.2 MeV. (6.3) 

The sensitivity to a 2 n is weak and we used 

02 ^ = 0.95 , 0.7 <Cp< 1.2 MeV. (6.4) 

In addition to these parameters we found for the Unna 
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(5a = —0.0013 ; ca = 1-6 fm ^ , 


(6.5) 


for all Cp and Wq considered. 


B. Variational energies 

Tables I and II show various components of the ^)(0 energy for the cases of no ANN 
potential and Cp = 0.7 MeV, Wo = 0.015 respectively. The cluster expansion for terms 
involving the A is converging well and it appears that it is not necessary to extrapolate 
beyond the four-body clusters for these terms. To get an accurate total energy of the ^)(0 
nucleus, it would be necessary to extrapolate the values of Vnnn as was done in PWP. 
However since we are mainly interested in Ba, we do not do that here and instead subtract 
an unextrapolated energy. 

One important result shown in Table II is that the expectation of is substantial 
and negative. This arises from the non-central correlations in the wave function. A purely 
Jastrow wave function {Uij = = UiA = flpA = 0) gives (Vanja) — 0-9 MeV; including Uij 

(with a = 0.94), Uijk, and UiA (with 0:277 = 1-0), but no Uja, results in = -4.2 MeV. 

Including UijA also but still keeping o = 0.94, 02,7 = 1.0, gives (I^atat) = -7.9 MeV, while 
lowering o to the optimal value of 0.87 and 02,7 to 0.95 reduces this to -5.0 MeV (this loss 
of binding is offset by changes in the expectation values of other parts of the Hamiltonian). 
These results can be understood as follows: by using the relation 

<T-Acr-B = A- B-|-i<T-(AxB), (6.6) 

the {Via, V 2 a} appearing in can be expressed in terms of the operators cr^ ■ fia cr 2 ■ r 2 A, 
(Ti ■ riA cr 2 ■ riA, cti ■ r 2 A (^2 ■ r 2 A, and cti ■ (T 2 and hence is a generalization of the tensor 
operator Su- However the expectation value of S '12 in a Jastrow wave function for a closed- 
shell nuclear system is zero, while the expectation value of 8^2 is non-zero. Hence the 
S '12 operator in U 12 completely changes {Va 12 )- Of course the UijA further enhances its 
contribution. 


15 



The energy that corresponds to the present calculation {Uijk acts first, use of only 
the hrst six operators in Argonne ^ 14 , and no extrapolation) is -101.0(9) MeV. We emphasize 
that, because of the above approximations, this energy is to be used only in comparison with 
the ^)(0 energies. The resulting i?A(^AO) 27.5(2.0) MeV for no Vaatat and 13.5(1.8) MeV 
for Cp = 0.7, ITo = 0.015, which are to be compared with the empirical value of 13.0(4) found 
in Sec. II. Thus even with realistic NN potentials and correlations, is very overbound if 
no Vavat is used. However a reasonable Vann results in a Ba consistent with the empirical 
value. 

To study the dependence of Ba on the strength of the Vann, we made a number of cal¬ 
culations with different values of Cp and Wq. In each case the NN quenching parameter was 
chosen according to Eq. (6.3). To minimize statistical errors, we made correlated difference 
calculations using the Cp = 0.7, Wq = 0.015 random walk. Table III presents the resulting 
changes, SBa, in Ba- All of these values of Ba and SBa are well £t by the formula 

Ba = 27.3 - 8.9 Cp + 11.2 C^ - 870. Wo; (6.7) 

the statistical error of Ba is ± 1.6 MeV. The quadratic dependence on Cp comes from 
UijA', the contribution of the dispersive term in UijA (and also in Uijk) is very small. When 
comparable calculations of other hypernuclei are available, Eq. (6.7) and the empirical value 
of Ba(^aO) = 13.0(4) can be used to uniquely determine the values of Cp and Wq (or to 
show that a different Hamiltonian is needed if no fit can be found). 

C. Densities and polarization of core 

Figure 2 shows point nucleon and A densities for the calculations of Tables I and H. The 
density of is also shown. For the case with no Vann, the nuclear correlation parameters 
were not changed from those used in The resulting density is, however, reduced 
near the origin and somewhat more peaked at r = 1.4 fm. This is presumably due to the 
repulsive Jan which pushes the nucleons away from the A which is strongly localized near 
the origin. 
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With Vann, the NN quenching parameter was signihcantly reduced. This results in a 
slightly more repulsive fc {a does not quench the central part of Vnn) and so the nuclear 
density is reduced for r ^ 2.2 fm. The nuclear kinetic and potential energies (see Tables I 
and II) for the no Vann case are separately larger in magnitude due to the higher density. 
It is probably accidental that the total nucleon energies for the two cases are so nearly the 
same: -90(2) MeV. This value is 11 MeV less than the corresponding binding energy of 
showing that the A significantly reduces the binding of the nucleons. 

The density prohle of the A for the two cases is also shown in Fig. 2, along with the density 
corresponding to the one-body part of T, i.e. |0A(r)p/47r. In both cases the Jastrow part 
of T, i.e. n/^(rjA), signihcantly increases the A density at the origin. This is presnmably 
because is small for r^A —^ 0 and hence pushes the A away from the high nuclear density 
around r = 1.4 fm. Because this density is larger for Vaattv = 0, the central A density is also 
larger in this case. 


VII. CONCLUSIONS 

In this stndy, we have extended the cluster Monte Carlo techniqne developed in PWP 
for to The cluster contributions that involve the A seem to converge well. It 
thus seems sufficient to include terms up to four-body clusters in the calcnlation. These 
calculations have been performed for a nnmber of sets of Wq and Cp which will be helpfnl in 
determining the parameters of the three-body interaction U^a by htting the Ba valnes of ^)(0 
and other hypernuclei. The present calculations show that the use of non-central NN, NNN, 
NA and NNA correlations completely change the expectation value of the three-baryon ANN 
interaction fonnd with central wave fnnctions, and thus have a strong effect on the choice 
of the parameters Wq and Cp. With snch correlations, reasonable valnes of Cp and Wo give 
the correct i?A(^AO)- ^Iso hnd that the A signihcantly changes the density prohle and 
energy of the 16 nucleons in ^)(0; the ANN potential is particularly signihcant in this regard. 
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FIGURES 


FIG. 1. Terms contributing to '^ann- 


FIG. 2. One-body densities for the nucleons and A in and for The short-dashed curve 
is the density corresponding to just (j)\. 
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TABLES 


TABLE I. Variational energies for Cp = Wq = 0. All values are in MeV. Numbers in paren¬ 
theses are statistical errors in the last digit. 


Clusters 


one-body 

two-body 

three-body 

four-body 

Total 

Kinetic Energy^ 


20.2(6) 

0.1(1) 

1.1(5) 

-0.7(3) 

20.7(9) 

AN Potential 

- e) 


-44.9(9) 

3.7(3) 

1.4(7) 

-47.3(12) 


vo{r)ePx 


-13.4(4) 

0.4(1) 

0.7(3) 

-12.2(5) 


jVaT^(TA ■ (TN 


0.34(3) 

-0.06(1) 

-0.05(2) 

0.22(3) 

A Energy 


20.2(6) 

-57.9(11) 

-2.3(7) 

1.4(11) 

-38.6(10) 

Nuclear Kinetic 


317.(2) 

269.(2) 

-19.(3) 

11.(3) 

556.(5) 

NN Potential 

Vij 


-737.(4) 

111.(2) 

4.(4) 

-623.(5) 

NNN Potential 

^ijk 



-59.8(8) 

35.8(8) 

-23.9(8) 

Nuclear Energy 


317.(2) 

-468.(2) 

32.(2) 

29.(3) 

-90.(2) 

Total Energy 


337.(2) 

-526.(3) 

30.(3) 

30.(1) 

-128.5(18) 


^Includes nucleon kinetic energy from AN correlations. 
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TABLE II. Variational Energies for Cp = 0.7 MeV, Wq = 0.015 MeV. All values are in MeV. 


Numbers in parentheses are statistical errors in the last digit. 


Clusters 


one-body 

two-body 

three-body 

four-body 

Total 

A Kinetic Energy'^ 


17.1(5) 

0.1(1) 

2.4(5) 

-1.0(4) 

18.6(6) 

AN Potential 

uo(r)(l - e) 


-39.5(7) 

-2.0(2) 

0.6(5) 

-40.9(8) 


vo{r)ePx 


-12.5(3) 

0.5(1) 

-0.2(4) 

-12.2(4) 


jVaT^CTA ■ ctn 


0.26(2) 

-0.06(1) 

0.00(2) 

0.20(2) 

ANN Potential 

'^ANN 



14.1(4) 

-0.1(1) 

14.0(4) 


T/27r 

^ANN 



-6.5(3) 

1.5(2) 

-5.0(3) 

A Energy 


17.1(5) 

-51.7(10) 

8.4(6) 

0.9(9) 

-25.3(7) 

Nuclear Kinetic 


309.(2) 

229.(2) 

-7.(2) 

-12.(2) 

520.(4) 

NN Potential 

Vij 


-682.(3) 

85.(2) 

10.(3) 

-587.(4) 

NNN Potential 

^ijk 



-50.2(6) 

27.9(6) 

-22.3(7) 

Nuclear Energy 


309.(2) 

-453.(2) 

28.(2) 

27.(2) 

-90.(2) 

Total Energy 


326.(2) 

-505.(2) 

36.(2) 

38.(2) 

-114.5(16) 


^Includes nucleon kinetic energy from A correlations. 
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TABLE III. Differences B\{Cp, Wo)-Bt^{Cp = 0.7, lEo = 0.015). The NN quenching parameter, 


a, is also given. 


Cp 

hEo 

a 

(5Ba MeV 

0.7 

0.01 

0.87 

+4.4 ± .3 

0.7 

0.017 

0.87 

-1.7+ .1 

0.7 

0.02 

0.87 

-4.4 + .3 

0.9 

0.01 

0.85 

+6.3 + .7 

0.9 

0.015 

0.85 

+1.9 + .5 

0.9 

0.02 

0.85 

-2.4 + .5 

1.0 

0.015 

0.84 

+2.8 + .8 

1.0 

0.02 

0.84 

-1.5 + .7 
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